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The primordial curvature perturbation £ may be generated by some curvaton field a, which is 
negligible during inflation and has more or less negligible interactions until it decays. In the current 
scenario, the curvaton starts to oscillate while its energy density p a is negligible. We explore the 
opposite scenario, in which p a drives a few e-folds of inflation before the oscillation begins. In this 
scenario for generating £ it is exceptionally easy to solve the n problem; one just has to make the 
curvaton a string axion, with anomaly-mediated susy breaking which may soon be tested at the 
LHC. The observed spectral index n can be obtained with a potential V oc <f> p for the first inflation; 
p = 1 or 2 is allowed by the current uncertainty in n but the improvement in accuracy promised 
by Planck may rule out p = 1. The predictions include (i) running n' ~ 0.0026 (0.0013) for p = 1 
(2) that will probably be observed, (ii) non-gaussianity parameter /nl ~ — 1 that may be observed, 
(iii) tensor fraction r is probably too small to ever observed. 

PACS numbers: 98.80Cq 



Introduction — The primordial curvature perturba- 
tion £ is already present a few Hubble times before cos- 
mological scales start to enter the horizon. At that stage, 
and on those scales, its Fourier components £(k) are time- 
independent and set the principle (or only ) initial con- 
dition the subsequent formation of large-scale structure 
in the Universe As a result, £(k) can be determined, 
and one of the main tasks of theoretical cosmology is to 
explore models of the early universe that can generate it. 

The generation of C(k) presumably starts at horizon 
exit during inflation {k = aH = a where a(t) is the scale 
factor of the universe) when the vacuum fluctuation of 
one or more scalar (or vector Q) fields becomes a clas- 
sical perturbation. According to the original scenario, £ 
is generated by the perturbation Sip of the inflaton field 
in a single-field slow-roll inflation model. In that case, 
£ is generated promptly at horizon exit, remaining con- 
stant thereafter. According to the curvaton scenario 
£ is instead generated by the perturbation 8a of a 'curva- 
ton' field, that has practically no effect during inflation 
and generates £ only when it's energy density becomes a 
significant fraction of the total. 

Instead of the curvaton scenario one can consider an 
inflaton-curvaton scenario where both 8<p and 8a con- 
tribute significantly to C 01 ■ O ne can a l so suppose that 
£ is generated during multi-field inflation, or by a 'mod- 
ulating' field that causes an effective mass or coupling to 
be inhomogeneous [l[. In this Letter we stay with the 



simpler curvaton scenario. 

Up till now, it has been assumed (Figure 1) that the 
curvaton starts to oscillate while its energy density p a is 
a negligible fraction of the total. Here we assume instead 
(Figure 2) that p a comes to dominate the total while it 
is still slowly varying, giving rise to a second era of infla- 
tion. For simplicity, we demand that these 'cosmological 
scales' are outside the horizon (k < aH) when the second 
inflation begins. 

A second inflation has been discussed many times be- 
fore, but always within a scenario where the first inflation 
generates at least a significant part of £, and/or some or 
all cosmological scales start out within the horizon. The 
second inflation is usually supposed to begin while the in- 
flaton of the first inflation is still oscillating. The contri- 
bution of the second inflation for this 'double inflation' is 
calculated for instance in [|| taking cosmological scales 
to be outside the horizon, and in [7| taking them to be 
partially inside. A second inflation starting during radi- 
ation domination might be called late inflation. Slow-roll 
late inflation was considered in [§], and fast-roll late in- 
flation in 0, but these authors ignore the effect on C of 
the second inflation. There is also thermal inflation [lfj 
which really has no effect on £. Our scenario is different 
from all of these. 

Duration of the second inflation — By 'cosmolog- 
ical' scales we mean those probed more or less directly 
by observation. They range from k = ciqHq where de- 
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FIG. 1: Original curvaton scenario. The inflaton <f> gives 
nearly constant energy density during inflation which af- 
ter reheating (shown as instantaneous) converts to radiation 
with p T oc a -4 . The curvaton density p a varies slowly with 
pa <C pi, until a starts to oscillate giving p„ a a" 3 . It is 
assumed that C, is constant after the curvaton decays, which 
is guaranteed if the universe is then radiation dominated as 
in the Figure. 
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FIG. 2: Inflating curvaton scenario. Here p a exceeds p r be- 
fore a oscillates, giving a few e-folds of inflation before the 
oscillation starts. 



notes the present, to around k = e 15 ao-ffo- According to 
a standard calculation [l| , requiring that they are outside 
the horizon when the second inflation starts corresponds 
to 



N 2 < 45 - ln(Kr 5 Mp/i7 2 )/2, 



(1) 



where N 2 = A (In a) is the e-folds of expansion during 
the second inflation, H 2 = a/ a is the Hubble parameter 
then and M P = (87rG)" 1/2 = 2 x 10 18 GeV. [A subscript 
1(2) will always denote the first (second) inflation.] 

The right hand side of Eq. takes p = 3M$H 2 to 
be constant during the second inflation and then oc a~ 4 
(radiation) till the observed matter-dominated era. Re- 
placing some of the radiation domination by matter dom- 
ination reduces it. The second term is positive because 
H% < Hi and we need Hi < 10~ 5 Afp or the tensor per- 
turbation with spectrum Vh = (8/Afp)(_ffi/27r) 2 would 
have been observed (see Eq. © below). 



CDM and baryon number cannot be created before 
the curvaton (or any other mechanism) creates £, since 
that would give an isocurvature perturbation excluded 
by observation Q. That requires something like p 2 ^ 4 > 
10 3 GeV corresponding to H 2 > KT 30 M P and N 2 < 16. 
We therefore require roughly N 2 < 45 to 16. 

Calculating the curvature perturbation — The 
curvature perturbation £ is described non-perturbatively 
through the SN formalism as in 11]. In this paper, we 
just work to first order in Sp, as in [3J. Then 



C(M) = -#(*) 



5p(k,t) 1 6p(k,t) 



(2) 



P (t) 3 P (t)+p(ty 

where p is the pressure, and Sp is defined on the slic- 
ing with uniform locally-defined scale factor a(x, t) (flat 
slicing). The second equality corresponds to the energy 
continuity equation p = —3H(p + p). 

Keeping only super-horizon scales, the e nerg y continu- 
ity equation is valid locally. As a result [12], C is con- 
stant during any era when p(x, t) is a unique function 
of p(x, t); that is the case for pure radiation (p = p/'i) 
or matter (p = 0). When cosmological scales start to 
enter the horizon, the temperature T is somewhat less 
than 1 MeV and we know that the Universe is practically 
pure radiation giving a time-independent £(k, t) that we 
are denoting simply by C(k). According to the curvaton 
scenario, C(k, t) does not vary between curvaton decay at 
t = td and T ~ MeV, which is guaranteed if the universe 
is completely radiation-dominated throughout that era. 

In any curvaton scenario, Q is generated while p = 
p a + p r and p = p a + p T , where p a = V{cr) + & 2 /2 and 
p a = —V(a) + a 2 1 2 are the curvaton contributions. For 
the original curvaton p v cannot be matter and is taken to 
be pure radiation, p r oc a~ 4 For the inflating curvaton, 
p r might be matter, p r oc a~ 3 allowing it to correspond 
to the oscillation of a field <j> responsible for the first in- 
flation. It might even decrease more slowly, say like a -1 
corresponding to a cosmic string network 13]. In these 
cases though, we demand that p r becomes radiation be- 
fore it is a significant fraction of p, so that it does not 
cause C to vary significantly. 

To facilitate an analytic calculation, one writes Eq. ([2]) 
as C(k,£) = /(i)Ca(k,t) where 3C CT = (8p a )/(p a +p a ) and 
/(<) = [p a +p a )/{p + p). There is supposed to be neg- 
ligible exchange of energy between the two components, 
so that Co- is constant if p a (x,t) is a unique function of 
/Ccr(x, t). For the original curvaton, £ CT becomes constant 
only after the oscillation begins, when p a ~ 0. We now 
argue that for the inflating curvaton, £ a will become con- 
stant soon after the second inflation begins at the epoch 
t = t 2 . We begin with the following equation, valid in 
the absence of perturbations [l[ : 



a + 3H& + V'{a) = 0. 



(3) 



Since sub-horizon modes of er(x, t) redshift away that 
quantity has negligible spatial gradient [ij . As a result it 
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satisfies Eq. Q for each x, with 3MpH 2 = p a + p T the 
locally defined quantity and t the proper time. Also, p a 
is a unique function of p a if and only if <r(x, t) is unique 
up to the choice of t — (attractor solution) . We expect 
that soon after the inflation begins it will be slow-roll or 
fast-roll (see below), and both of these make a a unique 
function of a giving indeed the required attractor solu- 
tion. 

If p T were completely negligible we could invoke a more 
general argument for the attractor but that might 
not apply because although the contribution of p T to H 
becomes small soon after the second inflation begins, its 
contribution to 



2e H = 2\H\/H Z =3 <j 



3 P > 



IP 



(4) 



may be dominant at least initially. While that is happen- 
ing, fit) ~ (d 2 )/(4jO r /3) <C 1. But just before a decays 
at td we have f(td) — Pal\Pa + (4/3)p r ] which will be 
very close to 1, In contrast, the oscillating curvaton can 
have /(id) <C 1. 

Keeping only super-horizon modes, f(td) = 1 gives 
C(x) = Sp a (x,t 2 )/3a 2 {t2). To first order in 5a, C(x) = 
V'5a(K)/3a 2 . At horizon exit during the first inflation 
5a is nearly gaussian with spectrum Hi/2tt. Allowing 
<r(x, t 2 ) to be a function g of its value at horizon exit, 
but taking both <r(x, t) and Hi to be time-independent 
while cosmological scales leave the horizon, we get the 
scale-independent spectrum 



>V2 „ 9 



' V>(t a )) Hx 



3 a%(t 2 ) 2tt 



(5) 



Observation gives 15] Vl /2 ~ 5 x 1CT 5 . 

Taking into account the time-dependence of H\ and a 
while cosmological scales leave the horizon one finds [ij 



i(k) - 1 = d\nV/d\nk = -2e H1 + 2r)i, 



(6) 



with the right hand side evaluated at horizon exit k = aH 
and 3H 2 n = d 2 V/da 2 . [Note that V = V(a, <j>, ■ ■ •) 
during the first inflation where </>,••• are the fields re- 
sponsible for that inflation.] Assuming a tensor fraction 
r = Vh/Vc. < 1CT 1 and \n'\ < 10" 1 where ri = dn/hik, 
observation gives n- 1 = -0.037 ± 0.014. 

Eq. © is a universal formula, applying whenever a 
field a different from the inflaton generates £. It can 
easily happen (as in our case, see below) that the last 
term is negligible. Then we need 2chi — 0.037, leading 
to three important consequences. 

1. To get the required eni we need [l6| a large change 
in the inflaton field <f> during the first inflation. To achieve 
that one usually takes <j) to have the canonical kinetic 
term with V(<t>, a) ~ V(<p) = A<$P . Then e H i ^ p/4A r (fc), 
where N(k) is the number of e-folds of the first infla- 
tion after the scale k leaves the horizon [jj. Defining 



Ni = N(a H ), we need Ni ~ 14p to get the required 
2chi — 0.037. While one is free to postulate any p the 
only choices of p with good justification are p — 1 ( cor- 
responding to monodromy |17j ) and p = 2 {corresponding 
to ' extranatural' [l8| inflation or N-flation [l9| . These 
give Ni ~ 14 and 28 respectively. 

2. A standard calculation [l| gives Ni ~ 60 — 
ln(10 _5 M P /iJi)/2 - N 2 , where the equality would be 
exact if p were constant during both inflations and cx 
a~ 4 otherwise until the observed matter-dominated era. 
Combining this with Eq. ([T]) gives 

15 + | ]n{Hi/H 2 ) < N t < 60 - 1 \n{lQ- h M v /H x ). (7) 

Taking into account the uncertainty in n this is compat- 
ible with N ~ 14p for p = 1 or 2, though p = 1 may 
be ruled out when Planck reduces the uncertainty. For 
the oscillating curvaton, where N 2 is absent, we would 
probably need p ~ 3 to 4. 

3. Since n(k) — 1 cx 1/N(k) the 'running' n' is given 
by n' = (1 — n)/Ni. This prediction holds also for the 
oscillating curvaton scenario (if |^i| e_ffi), and for the 
inflaton scenario within some simple slow-roll models and 
it makes n' big enough to observe in the future. For 
the oscillating curvaton, and the inflaton scenario, one 
expects roughly Ni ~ 60 corresponding to n' ~ 0.0007. 
For the inflating curvaton, we have for V((j)) cx (j? and 
taking account of the uncertainty in n 



1-n 
0.037 



0.0026 
P 



(8) 



with p = 1 or 2. This precise prediction for n' will be 
probably be tested in the future ( see for instance 20] ) . 

We require the contribution of the first inflation to be 
negligible, s = Vq^/Vq 1. Assuming canonical kinetic 
terms for the inflaton(s) [l|, 



V. 



1/2 



> 



1 



Hi 



V2TJH 27rAf P 



(9) 



and r < l&stui implying H\ < 1.1 x 10 _4 (rs) 1 / 2 Mp. 
(The equalities hold for a single inflaton.) Observation 
15] gives r < 10 _1 . Even if the second term of Eq. ([6]) 
contributes significantly, it is unlikely to accurately can- 
cel the first term which means that we need chi ^ 0.02, 
giving r < 0.3s and H x < 6 x 10" 5 s 1/2 M P . The tensor 
fraction will not be observed by Planck 1211 if s < 10 _1 , 
and it will probably never be observed [22| if s < 10 -3 . 

To calculate the non-gaussianity parameter /nl we ex- 
pand 5p a to second order in 5a giving 



C(x) = V'5a/3a 2 + (3/5)f NL (V'5a/3a) 2 
/nl = {5a 2 V"/V' 2 ) (1 ' "'"-'^ 



9"l9' 2 ) 



(10) 
(11) 



Observation [15J requires — 10 < /nl < 74 which means 
that the second term of Eq. (jTU)) gives a negligible con- 
tribution [23[ to Vq. We need |/nl| ft 1 if /nl is ever 
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to be detected. With such a value it will indeed be a 
ood approximation to ignore the non-gaussianity of 8a 
24 | . But our first-order treatment of 8p a is reliable only 
for |/nl| ^ 1; for |/nl| ~ 1 ° ne should go to second 
order, or use the SN formalism as in Q. 

The usual way of achieving inflation would be through 
the slow roll approximation a ~ — V jZHi- That will 
typically make N 2 too large but let us anyway see what it 
implies. Differentiating it requires e# 2 <C 1 and 1 773 j <C 1 
where 772 = V" /3H 2 ■ Using Eq. (j4]) the former condition 



requires £2 
get 



M*(V'/V) 2 /2 < 1 and p T < p. Then we 



,1/2 _ 



(2e 2 )- 1 / 2 (H 1 /2 7 :)(g'/M P ) 



(3/5)/ NL = m(i + g"/g' 2 ) 



(12) 
(13) 



We would need \g"\ ^> g' 2 to get a detectable /nl- 

Most discussions of the oscillating curvaton take 
V(a) ~ m 2 r a 2 /2. This choice is impossible for the in- 
flating curvaton scenario if the first inflation has in- 
flaton(s) with canonical kinetic terms. Indeed, using 
Eqs. ([9]) and (TTBI the curvaton contribution V^ a is given 
by V^/Vq = 2chi iV 2 (ct 2 /o - ! ) 2 < 1 which means that 
VQ a cannot dominate. 

The curvaton a string axion — Any scheme for gen- 
erating £ from some field x encounters the r\ problem; 
that a generic supergravity theory gives pBj in the early 
universe \r}™\ ^ 1, where rf£ = (d 2 V/dx 2 )/3H 2 1 and 
3MpH 2 r = p aI excludes any radiation contribution. It is 
a problem for two reasons; (i) we need |?7 X | < 10 -2 while 
cosmological scales leave the horizon to keep \n— 1| small 
enough, (ii) unless |c?V"/c?xl is exceptionally small we gen- 
erally need \rj x \ <C 1 at all times or x wm be quickly 
driven to a minimum of V. 

For the inflaton scenario, the r\ problem exists only 
during inflation when H ~ const; it is often ignored and 
could be regarded as a fine-tuning requirement on the 
parameters of the supergravity inflaton potential. For the 
oscillating curvaton scenario the r\ problem may be more 
severe because it may exist for a long time after inflation, 
with H strongly decreasing. For the inflating curvaton 
it is definitely more severe because it exists during both 
inflations with very different values for H. 

To avoid the rj problem for the curvaton one can take 
it be a pNGB with the potential 



V(a) = \v 



1 + cos ( 

O~0 



(14) 



ractically independent of other field values. It is known 
that do <C Mp gives the oscillating curvaton, but it is 
unclear how to motivate such a value. Choosing instead 
do Mp would give a second inflation with 1 772 1 1, 
but N2 would typically be too big and this choice anyhow 
seems impossible within string theory (2(| . (If the latter 
difficulty is ignored we can replace a in Eq. (|14[) by <f> to 
arrive at 'Natural Inflation' llf.) 



What we need for the inflating curvaton is Co ~ Mp. 
Then, in the regime a <C Co we have V ~ Vq — m 2 a 2 /2 
giving 1 772 1 - — ' 1 and just a few e-folds of inflation. Setting 
3M 2 iJ 2 2 = V gives H^/m 2 ~ (2/3ir*)(o%/M$). The 
required value Co ~ Mp is achieved if a is a string axion 
with gravity- or anomaly mediated susy breaking. [26j |. 
Then m is of order the gravitino mass m g 26]. The 
curvaton and gravitino have to decay before they can 
upset BBN which requires [53] m > 10 4 GeV - 10 _ 14 M P . 
This corresponds to anomaly mediation, which gives m ~ 
10 4 to 10 5 GeV and (like any version of susy) may soon 
be tested at the LHC. 

The contribution of 771 to n — 1 will be negligible if 
Hx/m > [3(1 - n)/2]" 1 / 2 = 4. With such a low m this 
is a mild requirement which we will take to be satisfied, 
so that Eq. © holds. 

Writing V' = —m 2 a, and taking H2 to be constant, 
Eq. ([3]) gives during the second inflation 



a ~ FH20-. 




m/H 2 . (15) 



The slow-roll regime is m <C H2 but we are interested 
in the 'fast-roll' [28| regime m > H2 corresponding to 
F > 1. (The approximation F ~ m/ H2 is adequate for 
m 2 /iff > 3.) In the fast-roll regime, Eq. (TTSl) is self- 
consistent if N2 3> 1 and <r < do. 

Since Eq. © is linear, g is a linear function and we get 



1/2 = l 



Hi 



FH2 J 2nai 6no~i 



(3/5)/ NL = -{FH 2 /mf 



-1. 



(16) 
(17) 



The result for /nl may be strongly modified by the cor- 
rection of second order in 8p a , but barring a strong can- 
cellation it seems that /nl may eventually be detectable. 

Since inflation ends at a ~ cto we have FN% ~ 
ln(Mp/o- 2 ). Since N 2 > 1, we have F < FN 2 , hence 
m/H 2 < FN 2 . Using Eq. $Vo$ with H x > H 2 > 10 4 GeV 
gives FN 2 < 24. Going the other way Eq. (flU)) with 
Hi - l(T 6 Afp implies FN 2 < 4. 

Eqs. (fTB|) and (fT7|) are roughly the same as those of 
the oscillating curvaton model with /(id) — 1- But the 
result for n' is different for the oscillating curvaton [![; 
there we might have n — 1 ~ 2rji with n' negligible, and 
even if Eq. (|SJ) holds we expect p ~ 4. 

One may worry about the assumption that a(t±) is 
near the top of the potential, given that the first inflation 
may be of long duration. For a given -Hi(i), the late- 
time probability distribution of a at the end of the first 
inflation can be calculated [2^] . Taking that distribution 
to apply and also taking Hi to be constant, we would 
need Hf > Vq to have a significant probability that a 
is near the top. This requires H\jH-x > MpjH\ > 10 5 
and (since H 2 > 10 4 GeV) H x > l(T 7 Afp. The former 
bound would probably make Eq. ([7]) incompatible with 
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p = 1. But to know whether the estimate Hf > Vb 
is realistic one would have to calculate the evolution of 
the probability distribution with the correct H\ (t) and a 
range of initial values of a. 

Conclusion — The hypothesis that the curvaton is a 
string axion leads to a simple early-universe scenario. 
The curvaton generates a few e-folds of inflation with 
H ~ 10 4 GeV, during which £ is created. The main in- 
flation takes place earlier, with a potential V{4>) oc (j) p , 
with p = 1 or 2 needed to reproduce the observed spectral 
index within current observational uncertainity. The hy- 
pothesis requires low-energy susy with anomaly-mediated 
susy breaking which may soon be tested at the LHC. It 
predicts that a tensor fraction r that is probably too small 
ever to observe, but a running n' that eventually be ob- 
served and will decide between the linear and quadratic 
potentials. A third prediction /nl ~ — 1 may also be 
testable, but the accuracy of the calculation needs to be 
improved. 
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